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Abstract 

In this paper, we address whether the growth of supermassive black-holes has kept pace with the process of 
galaxy assembly. For this purpose, we first searched the Hubble Ultra Deep Field (HUDF) for "tadpole galaxies" , 
which have a knot at one end plus an extended tail. They appear dynamically unrelaxed — presumably early- 
stage mergers — and make up ~6% of the field galaxy population. Their redshift distribution follows that of field 
galaxies, indicating that — if tadpole galaxies are indeed dynamically young — the process of galaxy assembly 
generally kept up with the reservoir of field galaxies as a function of epoch. 

Next, we present a search for HUDF objects with point-source components that are optically variable (at the 
> 3. Oct level) on timescales of weeks-months. Among 4644 objects to i^ s ~28.0 mag (IOct), 45 have variable point- 
like components, which are likely weak AGN. About ~1% of all field objects show variability for 0.1 < z < 4.5, and 
their redshift distribution is similar to that of field galaxies. Hence supermassive black-hole growth in weak AGN 
likely also kept up with the process of galaxy assembly. However, the faint AGN sample has almost no overlap 
with the tadpole sample, which was predicted by recent hydrodynamical numerical simulations. This suggests 
that tadpole galaxies are early-stage mergers, which likely preceded the "turn-on" of the AGN component and the 
onset of visible point-source variability by > 1 Gyr. 
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1. Introduction 

From the WMAP polarization results (Kogut 
et al. 2003), population III stars likely existed at 
z~20. These massive stars ( > 250M Q ) are ex- 
pected to produce a large population of black 
holes (BH; M 6fe >150M Q ; Madau & Rees 2001). 
Since there is now good dynamical evidence for 



the existence of supermassive (Mbh— 1O 6 -1O 9 M0) 
black holes (SMBH's) in the centers of galaxies 
at z~0 (Kormendy & Richstone 1995; Magorrian, 
Tremaine, & Richstone 1998; Kormendy & Geb- 
hardt 2001), it is important to understand how the 
SMBH's seen at z~0 have grown from lower mass 
BH's at z~20. A comprehensive review of SMBH's 
is given by Ferrarese & Ford (2004). One sugges- 
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tion is that they "grow" through repeated merg- 
ers of galaxies which contain less massive BH's, so 
the byproduct is a larger single galaxy with a more 
massive BH in its center. The growth of this BH 
may then be observed via its AGN activity. If this 
scenario is valid, there may be an observable link 
between galaxy mergers and increased AGN activ- 
ity (Silk & Rees 1998). Therefore, studying this 
link as a function of redshift could give insight into 
the growth of SMBH's and its relation to the pro- 
cess of galaxy assembly. 

Recent numerical simulations addressed some 
long-standing issues in the dissipational collapse 
scenario by including previously-neglected ener- 
getic feedback from central SMBH's during the 
merging events (e.g., Robertson et al. 2005). They 
emphasize the relationship between the central BH 
mass and the stellar velocity dispersion, which con- 
firms the link between the growth of BH's and 
their host galaxies (di Matteo et al. 2005; Springcl 
et al. 2005ab). The present study provides obser- 
vational support for these models at cosmological 
redshifts. 



2. The Hubble Ultra Deep Field data 

The Hubble Ultra Deep Field (HUDF; Beckwith 
et al. 2005) is the deepest optical image of a slice 
of the Universe ever observed. It consists of 400 
orbits with the HST Advanced Camera for Sur- 
veys (ACS) observed over a period of four months 
in four optical bands (BVi'z'). These are supple- 
mented in the J-ff-bands with the Near-Infrared 
Camera and Multi-Object Spectrograph (NIC- 
MOS; Bouwens et al. 2004). The HUDF reaches 
~1.0 mag deeper in B and V and ~1.5 mag deeper 
in i' AB and z' AB than the equivalent filters in the 
Hubble Deep Field (HDF, Williams et al. 1996). 

A large number of galaxies in the HUDF ap- 
pear dynamically unrelaxed, which suggests they 
must play an important role in the overall pic- 
ture of galaxy formation and evolution. In partic- 
ular, we notice many galaxies with a knot-plus- 
tail morphology, which constitute a well-defined 
subset of the irregular and peculiar objects in the 
HUDF that is uniquely measurable. According to 



di Matteo et al. (2005), this morphology appears 
to represent an early stage in the merging of two 
nearly-equal mass galaxies. They are mostly linear 
structures, resembling the "chain" galaxies first re- 
ported by Cowie, Hu, & Songaila (1995). When 
more than two clumps come together, these ob- 
jects may be more akin to the luminous diffuse 
objects and clump clusters (Conselice et al. 2004; 
Elmcgrccn, Elmegreen, & Sheets 2004; Elmcgreen, 
Elmcgrccn, & Hirst 2004), or other types of ir- 
regular objects (Driver et al. 1995; van den Bergh 
2002). 

Since the HUDF data was observed over a period 
of four months, it also provides a unique opportu- 
nity to search for variability in all types of objects 
to very faint flux levels, such as faint stars, distant 
supernovae (SNe), and weak active galactic nuclei 
(AGN). From all objects detected in the HUDF, 
we therefore selected the subset of tadpole galaxies 
and variable objects, and analyzed their properties 
in the i^ B -band, where the HUDF images are deep- 
est and have the best temporal spacing over four 
months. Yan & Windhorst (2004b) discuss how the 
i^ B -selection result in a small bias against objects 
at z>5.5 in the high redshift tail of the redshift 
distribution. However, tadpole galaxies at z~5.5 do 
exist (e.g., Rhoads et al. 2005). Since most HUDF 
objects have z > 1.5, the i' AB -ba,nd images sample 
the rest-frame UV, where AGN are known to show 
more variability (Paltani & Courvoisier 1994). 

To address whether supermassive black-hole 
growth kept pace with galaxy assembly, we will 
present in this paper the redshift distribution of 
both tadpole galaxies and weak variable AGN in 
the HUDF, and compare these with the redshift 
distribution of the general field galaxy population. 



3. Tadpoles as Proxy to Galaxy Assembly 

The steps to select galaxies with the character- 
istic "tadpole" shape are described in Straughn 
et al. (2006). In short, objects of interest have a 
bright "knot" at one end and an extended "tail" 
at the other. Two different source catalogs were 
made to i' AB —28.0 mag using SExtractor (Bertin 
& Arnouts 1996): a highly deblended catalog con- 
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taining many point-like sources, including the 
knots of potential tadpole galaxies, and a low- 
deblending catalog containing extended sources, 
including the tadpole's tails. 

First, the knots of the tadpole galaxies were 
selected by setting an axis-ratio limit. "Knots" 
were defined from the highly deblended catalog 
with an axis ratio rounder than some critical value 
(6/a>0.70). "Tails" are elongated objects selected 
from the low-dcblcnding catalog with 6/a<0.43. 
Tadpoles were defined when a knot was within a 
certain distance of the geometrical center of a tail, 
namely <4a (in semi- major axis units of the tail). 
We also required that the knot be >0.1a from the 
tail's geometrical center, since we are searching for 
asymmetric objects, and want to eliminate upfront 
as many of the true edge-on mid-type spiral disks 
as possible. The tadpole candidates also must have 
the knot near one end of the tail, hence we selected 
only those tails and knots with a relative position 
angle A6><20°, as measured from the semi-major 
axis of the tail. This prevented including knots 
and tails that appear close together on the image, 
but are not physically part of the same galaxy. 

Our final sample contains 165 tadpole galaxies, 
a subset of which is shown in Fig. 1. These were 
selected from 2712 objects in the low-deblending 
HUDF catalog to i' AB =28.0 mag. Less than 10% 
of the selected tadpoles appear as normal edge-on 
disk galaxies. Fig. 2 shows a significant overabun- 
dance of knots near the end of the elongated diffuse 
structures (A0<1O°) as compared to randomly 
distributed knots. Hence, the majority of tadpoles 
are not just chance alignments of unrelated knots. 
Instead, we believe they are mostly linear struc- 
tures which are undergoing recent mergers. Their 
redshift distribution is shown in Fig. 3a-3b. 



4. Faint Variable Objects as Proxy to 
SMBH Growth 

Our HUDF variable object study is described in 
Cohen et al. (2006). Individual cosmic-ray (CR) 
clipped images and weight maps were used with 
multidrizzle (Koekemoer et al. 2002) to create 
four sub-stacks of approximately equal expo- 



sure times that cover 0.4-3.5 months timescales. 
These used the same cosmic-ray maps and weight 
maps as the full-depth HUDF mosaics. All four 
epochs were drizzled onto the same output pixel 
scale (O'.'030/pixel) and WCS frame as the origi- 
nal HUDF. Since we are searching for any signs 
of variability, we used a liberal amount of object 
dcblcnding in the SExtractor catalogs, which 
used a l.Ocr detection threshold and a minimum of 
15 connected pixels (i.e., approximately the PSF 
area) above sky. This allows pieces of merging 
galaxies to be measured separately, to increase the 
chance of finding variable events in point-source 
components. Since each of the four epoch stacks 
have half the S/N-ratio of the full HUDF, the sam- 
ple studied for variability contains 4644 objects to 
i^ B <28.0 mag (>10<r). 

The ACS/WFC PSF varies strongly with loca- 
tion on the CCD detectors, and with time due to 
orbital "breathing" of the HST Optical Telescope 
Assembly. Hence, we cannot use small PSF-sized 
apertures to search for nuclear variability, as could 
be done by Sarajedini et al. (2003a) for the much 
larger WFPC2 pixel-size and the on-axis location 
of the WFPC2 camera. Instead, we had to use total 
magnitudes of the highly deblended ACS objects. 
Even though our total flux apertures may encom- 
pass the whole galaxy, any variability must come 
from a region less than the 0'.'084 PSF in size, due to 
the light-travel time across the variability region. 

The four epoch catalogs were compared to each 
other, resulting in six diagrams similar to Fig. 4a, 
which show the change in measured total magni- 
tudes in matched apertures as a function of the 
full-depth HUDF flux. The flux-error distribution 
was determined iteratively for each pair of observa- 
tions, such that 68.3% of the points lie within the 
boundaries of the upper and lower l.Ocr lines that 
represent the Gaussian error distribution (Fig. 4a). 
In order to demonstrate the Gaussian nature of 
this error distribution at all flux levels, the Amag- 
data were divided by the l.Ocr model line, and his- 
tograms were computed for the resulting normal- 
ized Amag data at various flux- levels in Fig. 5. 
These histograms are well fit by normalized Gaus- 
sians with cr~ 1.0. The HUDF noise distribution is 
not perfectly Gaussian, but with 288 independent 



3 



exposures in the i^ B -band, the error distribution 
is as close to Gaussian as seen in any astronomical 
CCD application. Once the ±1.0ct lines were de- 
termined, we find all objects that are at least 3. Oct 
outliers. Most outliers in Fig. 5 at Amag > 3. Oct are 
due to object variability, after pruning large ob- 
jects without visible point sources which suffered 
from SExtractor dcblending errors. In Fig. 4a, we 
show the ±1(7, ±3(7, and ±5ct lines, along with the 
actual data. The choice of 3. Oct implies that we 
should expect 0.27% random contaminants. 

In total, wc find 45 out of 4644 objects that show 
the signatures of AGN variability. These are vari- 
able at the > 3. Oct level, have a compact region in- 
dicative of a point source, and are devoid of visible 
image defects or object splitting issues. Less than 
one of these 45 is expected to be a random contam- 
inant. In total, 577 candidates were rejected due 
to crowding or splitting issues, or due to the lack 
of a visible point source. Fig. 4b shows the num- 
ber of ct by which each object varied for each of 
the 6 possible epoch-pairs. The colored symbols are 
for the 45 "best" candidates. Another 57 objects 
were found that are "potentially" variable can- 
didates. The four-epoch light-curves for these 45 
variable candidates are shown in Fig. 6. Of these, 
49% were discovered from a single epoch-pair (usu- 
ally indicative of a global rise or decline as a func- 
tion of time in the light-curve), 43% in two epoch- 
pairs, and only 5% (2 objects) in 3 epoch-pairs. 
Further details are given in Cohen et al. (2006). In 
summary, the variability fraction on a timescale of 
few months (rest-frame timescale few weeks to a 
month) is at least 1% of all HUDF field galaxies. 

Since the HUDF is in the Chandra Deep Field- 
South (CDF-S, Rosati et al. 2002), there exists 
deep X-ray data. Within the HUDF, there are 16 
Chandra sources (Koekemoer et al. 2004,2006), 
and we detect four of these as variable in the 
optical. One of these is a mid-type spiral with 
i' AB ~2\.2A mag, that belongs to a small group of 
interacting galaxies. Two others are optical point 
sources with i^ B ~21.12 mag and ~24.79 mag, 
showing little or no visible host galaxy. Both have 
measured spectroscopic AGN emission-line red- 
shifts at z~3 (Pirzkal et al. 2004). The detection 
of 25% of the Chandra sources as optically vari- 
able in the HUDF data shows that the variability 



method employed here is a reliable way of finding 
the AGN that are not heavily obscured. 

The faint object variability in the HUDF is most 
likely due to weak AGN, given the timescalcs and 
distances involved. Strolger & Riess (2005) found 
only one moderate redshift SN in the HUDF, so 
SNe cannot be a significant source of contamina- 
tion in our sample. Several other possible source 
of incompleteness in the variability study must be 
addressed. Non-variable AGN, or AGN that only 
vary on timescalcs much longer than 4 months, 
or optically obscured AGN would not have been 
detected with our UV-optical variability method. 
Sarajedini et al. (2003ab) had two HDF epochs 5- 
7 years apart, and found 2% of the HDF to be vari- 
able. It is thus possible that our sampled times- 
scale shorter than 4 months missed a factor > 2 of 
all AGN — the ones variable on longer time-scales. 



5. The Redshift Distribution of Tadpole 
Galaxies and Faint Variable Objects 

We calculate photometric redshifts of all HUDF 
galaxies to i' AB =28.0 mag ( > 20ct) from the BVi'z' 
photometry using HyperZ (Bolzonella et al. 2000), 
plus NICMOS JH (Thompson ct al. 2005) and 
VLT ISAAC if -band images where available. 
When compared to published spectroscopic red- 
shifts for 70 CDF-S objects (Le Fevre et al. 2005), 
our photometric redshifts have an rms scatter of 
0.15 in S=(z phot -z spec ) I (l+z spec ) if all 70 objects 
are included, and 0.10 in 5 when we reject a few of 
the most obvious outliers. 

The redshift distribution of all HUDF galaxies 
(solid line in Fig. 3a) is as expected, with the pri- 
mary peak at z < 1.0 and a generally declining tail 
at z~4-5. These trends were also seen in the HDF 
field galaxies (Driver et al. 1998). A deficit of ob- 
jects is apparent at z~l-2 due to the lack of UV 
spectral features crossing the BViz(+JH) filters. 
Unlike the HDF, this deficit occurs because the 
HUDF does not yet have deep enough F300W or 
/7-band data. The resulting redshift bias, however, 
is the same for both tadpoles, variable objects and 
the field galaxy population, and so divides out in 
the subsequent discussion. Within the statistical 
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uncertainties in Fig. 3a, the shape of the tadpole 
galaxy redshift distribution follows that of the field 
galaxies quite closely. This suggests that if tadpole 
galaxies are indeed dynamically young objects re- 
lated to early-stage mergers, they seem to occur in 
the same proportion to the field galaxy population 
at all redshifts. Tadpole galaxies may therefore be 
good tracers of the galaxy assembly process. The 
ratio of the two redshift distributions N(z) and the 
resulting percentage of tadpole galaxies is plotted 
in Fig. 3b. Overall, the percentage of tadpole galax- 
ies is roughly constant at ~6% with redshift to 
within the statistical errors for the redshifts sam- 
pled (0.1<z<4.5). 

In Fig. 7a, we show the photometric redshift 
distribution for all HUDF objects with i' AB < 28.0 
mag, and for our best 45 variable candidates. Their 
redshift distribution follows that of the field galax- 
ies in general, i.e., there is no redshift where faint 
object variability was most prevalent. We plot in 
Fig. 7b the ratio of the N(z) for variable objects to 
that of field galaxies, and show that the weak vari- 
able AGN fraction is roughly constant at approxi- 
mately 1% over all redshifts probed in this study. 



6. Discussion and Conclusions 

The fact that about 6% of all field galaxies 
are seen in the tadpole stage is an important 
constraint to hierarchical simulations. Springcl 
et al. (2005ab) predict a tadpole-like stage ~0.7- 
1.5 Gyr after a major merger begins, suggesting 
that the tadpole morphology represents an early- 
merger stage of two galaxies of roughly compara- 
ble mass. If this 6% indicates the fraction of time 
that an average galaxy in the HUDF spends in an 
early-merger stage during its lifetime, then every 
galaxy may be seen temporarily in a tadpole stage 
for ~0.8 Gyr of its lifetime, and may have under- 
gone ~10-30 mergers during its lifetime (Straughn 
et al. 2006). More complex mergers involving mul- 
tiple components may lead to irregular/peculiar 
and train-wreck type objects, and the luminous 
diffuse objects or clump-clusters, which dominate 
the galaxy counts at faint magnitudes (Driver 
et al. 1998). Given that tadpoles only trace a cer- 



tain type and stage of merging galaxies, the above 
statistics are a lower limit on the number of all 
mergers. 

The question arises if tadpole galaxies and ob- 
jects with point-sources that show signs of variabil- 
ity are drawn from the same population. Among 
our 165 tadpole galaxies, none coincide with the 
sample of 45 variable objects or with the CDF-S 
X-ray sources S (Alexander et al. 2005). At most 
one or two of the variable candidates resemble the 
tadpole galaxies of Straughn et al. (2006). 

A factor of three of all AGN may have been 
missed, since their UV-optical flux was obscured 
by a dust-torus. In the AGN unification picture, 
AGN cones are two-sided and their axes are ran- 
domly distributed in the sky, so that an average 
cone opening-angle of w implies that a fraction 1— 
sin(w) of all AGN will point in our direction. If 
0t>~45° (e.g., Barthel 1989), then every optically 
detected AGN (QSO) represents 3 -4 other bulge- 
dominated galaxies, whose AGN reflection cone 
didn't shine in our direction. Hence, their AGN 
may remain obscured by the dust-torus. Such ob- 
jects could be visible to Chandra in X-rays or to 
Spitzer in the mid-IR, although the available Chan- 
dra and Spitzer data are not deep enough to detect 
all HUDF objects to AB~28 mag. 

Together with the factor of > 2 incompleteness 
in the HUDF variability sample due to the limited 
time-baseline sampled thus far, the actual fraction 
of weak AGN present in these dynamically young 
galaxies may thus be a factor of >6-8 x larger than 
the 1% variable AGN fraction that we found in the 
HUDF. Hence, perhaps as many as >6-8% of all 
field galaxies may host weak AGN, only ~1% of 
which we found here, and another > 1% could have 
been found if longer time-baseline had been avail- 
able. Another factor of 3-4 of AGN are likely miss- 
ing because they are optically obscured, The next 
generation of X-ray and IR telescopes (Windhorst 
et al. 2006ab) and longer optical time-baselines are 
needed to detect all weak AGN in the HUDF. 

Recent state-of-the-art hydrodynamical models 
(diMatteoet al. 2005; Springelet al. 2005ab; Hop- 
kins et al. 2005) suggest that during (major) merg- 
ers, the BH accretion rate peaks considerably after 
the merger started, and after the star-formation 
rate (SFR) has peaked. Their models suggest that, 
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for massive galaxies, a tadpole stage is seen typi- 
cally about 0.7 Gyr after the merger started, but 
~1 Gyr before the SMBH accretes most of its mass, 
which is when the galaxy displays strong visible 
AGN activity. Since the lifetimes of QSO's and 
radio-galaxies are known to be < (fewx 10 7 )-10 8 
years (Martini 2004; Grazian et al. 2004, Jakob- 
sen et al. 2003), these models thus imply that the 
AGN stage is expected to occur considerably (i.e., 
> 1-1.5 Gyr) after the early- merger event during 
which the galaxy is seen in the tadpole stage. 

The observed lack of overlap between the HUDF 
tadpole sample and the weak variable AGN sam- 
ple thus provides observational support for this 
prediction. Hopkins et al. (2005) have quantified 
the timescales that quasars will be visible during 
merging events, noting that for a large fraction of 
the accretion time, the quasar is heavily obscured. 
In particular, their simulations show that during 
an early merging phase — our observed tadpole 
phase — the intrinsic quasar luminosity peaks, but 
is completely optically obscured. Only after feed- 
back from the central quasar clears out the gas, 
will the object become visible as an AGN. This 
should be observable by Spitzer in the mid-IR as a 
correspondingly larger fraction of IR-selected ob- 
scured faint QSO's. To study the relation between 
galaxy assembly and SMBH growth in detail, we 
need deeper surveys at longer wavelengths with 
the James Webb Space Telescope (JWST; Wind- 
horst et al. 2006a). The JWST photometric and 
PSF stability are crucial for this, since many of our 
HUDF objects show significant variability of less 
than a few percent in total flux. 
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Fig. 1. i^ B -band mosaic of a subset of the HUDF tadpole galaxy sample. Stamps are 3 arcsec on a side. The vast majority 
of our tadpole sample contains the distinctive knot-plus-tail morphology. 
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Fig. 2. Distribution of position angle 8, measured from the major axis of the diffuse component, for all off-centered knots 
found within r<4a ( <2") from the center of an elongated diffuse object in the HUDF. There is a clear excess of knots near 
|A0|~O°, confirming the linear structure of most tadpoles. 




7 



500 
400 
300 
200 
100 




All UDF Galaxies 
Tadpole Galaxies 



h-i" 

I | 



2 3 4 

Photometric redshift 




2 3 4 

Photometric Redshift 



Fig. 3a (LEFT). Photometric redshift distribution of galaxies in the HUDF. The solid black histogram shows the redshift 
distribution of all HUDF field galaxies to i' AB =28.0 mag, while the dashed red histogram shows the redshift distribution 
of the tadpole galaxies, multiplied by 16 for best comparison. Fig. 3b (RIGHT). Percentage of total galaxies that are 
tadpoles vs. photometric redshift. Within the statistical errors, ~6% of all galaxies are seen as tadpoles at all redshifts. 
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Fig. 4a (LEFT). Magnitude difference between two HUDF epochs of all objects vs. i^ fl -band flux from matched total 
apertures. The ±lcr, ±3<r, ±5cr lines are shown. Blue points show the |Amag| <±1.0<r points used to normalize the error 
distribution. Large red points show the "best" 45 variable candidates from all six possible epoch combinations, many of 
which were seen at >3.0cr in two or more epoch combinations. Fig. 4b (RIGHT). Number of a that each object varies 
for all six possible epoch combinations. Colored symbols indicate the "best" sample of 45 variable candidates from Fig. 4a, 
that are unaffected by local image deblending issues or weight map structures. Each object appears six times in this plot. 
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Fig. 5. Gaussian nature of the HUDF total-flux error dis- 
tribution at all flux levels. The Amag data from Fig. 4a 
were divided by the best- fit model l.Ocr lines. Histograms 
for the indicated magnitude ranges are well fit by normal- 
ized Gaussians (parabolas in log space) with <r~1.0. The 
almost indistinguishable dashed and solid lines are for the 
best-fit a (indicated in the legends) and for assumed a=l 
Gaussians, respectively. Hence, Fig. 5 shows that 3. Oct really 
means 3.0±0.1c. All objects with Amag >3.0<r not affected 
by deblending errors are variable candidates. 
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Fig. 6. Light curves of the 45 best candidates with signs 
of optical point-source variability. The change in measured 
total flux (average— individual epoch) is plotted vertically, 
and the number of days since the first epoch is plotted 
horizontally. 
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Fig. 7 a (LEFT). Photometric redshift distribution of all HUDF field galaxies to i' AB <28.0 mag (solid line), and for the 
"best" variable candidates (red dashed line) multiplied by 60 X for best comparison. The redshift distribution of the variable 
objects follows that of field galaxies in general. Fig. 7b (RIGHT). Percentage of HUDF objects to i' AB <28.0 AB-mag 
showing variable point sources as a function of redshift. Within the statistical uncertainties, about 1% of all HUDF galaxies 
show point source variability over the redshift range surveyed (0 <z < 5). Hence, SMBH growth as traced by the weak AGN 
fraction keeps pace with galaxy assembly as traced by the tadpoles in Fig. 3b. 
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